Phase diagram of the dissipative quantum particle in a box 



J. Sabio 1 , L. Borda 2 - 3 , F. Guinea 1 , and F. Sols. 4 

1 Instituto de Ciencia de Materiales de Madrid (CSIC), Cantoblanco, 28049 Madrid. Spain 
2 Research Group "Theory of Condensed Matter" of the Hungarian 
Academy of Sciences, TU Budapest, Budapest, H-1521, Hungary 
3 Physikalisches Institut, Universitdt Bonn, Nussallee 12, D-53115 Bonn, Germany and 
4 Departamento de Fisica de Materiales. Universidad Complutense de Madrid. 28040 Madrid. Spain. 

(Dated: January 16, 2009) 

We analyze the phase diagram of a quantum particle confined to a finite chain, subject to a dissi- 
pative environment described by an Ohmic spectral function. Analytical and numerical techniques 
are employed to explore both the perturbative and non-perturbative regime of the model. For small 
dissipation the coupling to the environment leads to a narrowing of the density distribution, and 
to a displacement towards the center of the array of accessible sites. For large values of the dissi- 
pation, we find a phase transition to a doubly degenerate phase which reflects the formation of an 
inhomogeneous effective potential within the array. 
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I. INTRODUCTION. 

The problem of a quantum particle interacting with an 
environment deserves special attention since it has impli- 
cations in fundamental areas as quantum measurement 
theory, quantum dissipation and quantum computation, 
among others. A quantum particle interacting with an 
environment consisting of a continuum of degrees of free- 
dom, the Caldeira-Leggett model,- is actually the sim- 
plest model which can be used to study the destruction 
of quantum coherence and the emergence of classical be- 
havior in the framework of quantum mechanics. 

It is generally believed that the Caldeira-Leggett model 
captures the essential features of the behavior of more 
complicated open quantum systems. Therefore, its study 
as a toy model can be justified even if its not connected to 
any particular experimental realization. However, some 
extensions of this model are known to be relevant in 
real systems: for instance when analyzing dephasing in 
a qubit 2 - or in a dissipative Josephson junction* 3 - An- 
other application recently pointed out is the study of the 
decoherence induced in mesoscopic systems by external 
gates^^ This can be seen as a consequence of Caldeira- 
Leggett model reproducing the long time dynamics of 
particles interacting with Ohmic environments^ 

Three variations of this model have been particularly 
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well studied: (i) The dissipative two level system,, 
which is probably the most analyzed model in the context 
of quantum computation, being the archetype of a qubit 
in the presence of a bath^ (ii) A particle moving in a 
periodic potential ) 13 ! 14 in a magnetic field^ and in both 
a periodic potential and a magnetic field i 16 ' 17 The case 
of a dissipative particle in a periodic potential is relevant 
to the study of defects in Luttinger liquids^ while the 
case with both a magnetic field and a periodic potential 
applies to a junction between more than two Luttinger 
liquids. 19 Finally, (iii) the dissipative free particl o 20 ' 21 of 
interest in the study of quantum Brownian motion. In 
the first two cases the system undergoes a phase tran- 



sition, for a critical value of the coupling to the bath, 
to a phase where the particle is localized. This kind of 
transition, which belongs to the general class of bound- 
ary quantum phase transitions^ 2 - has been studied in the 
literature with many different approaches including path 
integral, renormalization group and variational Ansatz. 

In the present work, we study the phase diagram of a 
particle confined to a finite tight-binding chain coupled 
to an Ohmic dissipative environment through its coordi- 
nate variable. This is the simplest intermediate instance 
between two of the limiting cases described above, the 
dissipative two level system and a particle in an infinite 
array. It represents a quantum particle interacting with 
an ohmic reservoir whose motion is restricted to a fi- 
nite region. As discussed in detail below, the inclusion 
of the hard wall boundary conditions introduces inhomo- 
geneities in the density distribution of the particle, and it 
yields a non-trivial phase diagram, with a quantum crit- 
ical point which can be characterized in detail by using 
numerical techniques. 

The paper is organized as follows: First we describe 
the model, and briefly review the main results of related 
models obtained in the past. Then, we discuss the main 
techniques employed to analyze the ground state of the 
system. The calculated phase diagram is discussed, as 
well as the way various observable quantities are affected 
by the dissipation. The last part of the paper summarizes 
the main results of the work. 



II. THE MODEL. 



As we have noted above, the model we address is that 
of a particle coupled to a dissipative bath which can hop 
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between M sites. The Hamiltonian reads: 
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Here i is the hopping between nearest neighbor sites, and 
A the coupling with which the bath couples to the posi- 
tion of the particle q = X)m( m — m o)ctn c m, where mo 
labels the center of the chain. The bath is characterized 
by a high energy cutoff lo c . The last term of the Hamil- 
tonian is a counter-term introduced in order to preserve 
the degeneracy between the energies of the different sites. 
We assume that the coupling leads to Ohmic dissipation, 
so that: 

J{uj) — 7rA 2 kS(uj — k) — 2ira\uj\ , for u< lu c 

k<cj c 

(2) 

where a = A 2 / (4tt) and J {to) is the spectral weight func- 
tion. 

This Hamiltonian can be transformed through a uni- 
tary operation into one in which all sites are treated on 
the same footing because the transformed bath couples 
to the inter-site hopping. In order to arrive at this form, 
which can be useful in the study of the phase diagram, 
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we use the transformation U 
Hamiltonian, yielding: 13,23 
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Hamiltonians |1| and contain two limiting cases of 
interest. The case A = is that of a confined particle 
decoupled from the bath. The Hamiltonian is readily 
diagonalized, and the ground state corresponds to a par- 
ticle delocalized with a density p m = jAj sin 2 ( ™ ± ) . 
On the other hand, the case t = corresponds to a parti- 
cle without kinetic energy. From Hamiltonian ((3|) we see 
that the ground state is M-fold degenerate in the sub- 
space of site states adiabatically dressed by a cloud of 
bosons, |m) <§5 e ~ mX ^k Vk < - b >*'~ bk " , \oy l n general, we will 
be mainly interested in generic values of the parameters t 
and A of the Hamiltonian. Here, the bath can be regarded 
as performing repeated measurements of the position of 
the particle, localizing it in the sites basis, as opposed to 
the kinetic term, which tends to delocalize it. As in the 
dissipative two level system and the dissipative particle 
in a periodic potential, the phase diagram is expected to 



reflect the two opposing tendencies through a quantum 
phase transition. 

In this work we will analyze the regime i « w c , where 
we expect that the low-energy properties will depend only 
on the dimensionless parameters t/u c and a. Notice that 
in the continuum limit, M » I, where the couplings 
satisfy a << 1, the model should reduce to a particle 
described by an effective mass in a dissipative environ- 
ment, which admits a complete analytical solution^ For 
larger couplings we expect, as mentioned, a phase tran- 
sition as seen in related models. To analyze this region 
we will concentrate in a small number of sites, ranging 
between 2 and 6, where both numerical and analytical 
calculations are easier to perform. 



III. THE CALCULATIONS. 

A. Preliminary remarks. 

The existence of a phase transition for a large enough 
strength coupling can be seen simply by normal ordering 
Hamiltonian ([3]). This operation, which can be regarded 
as a resummation of an infinite series of tadpole diagrams 
for the interacting vertex^ gives a renormalized hopping 

tren = teT ~ * . Performing the summation we can 
write from this expression a flow equation for the hopping 
parameter in terms of the cutoff of the bath, giving the 
celebrated result^ 
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From this equation it follows that a phase transition ex- 
ists for a — 1 from a regime in which t ren is finite to one in 
which is effectively renormalized to zero, thus suppress- 
ing quantum fluctuations. In general, however, this is not 
the whole story and more precise calculations are needed 
in order to get the right critical line, which can have a 
dependence on the hopping parameter. This is the case, 
as is well known, for the dissipative two level system, 
where the transition can be shown to be of the Kosterlitz 
- Thouless (KT) kind 2 ^ when higher corrections to the 
flow equations are computedi&^i^ 3 - On the contrary, for 
the dissipative particle in a periodic potential there are 
no higher order corrections to the flow equations and the 
transition line in the a — t phase diagram can be shown 
to be vertical^ 

A similar analysis can be tried for the dissipative con- 
fined particle, but the computation of higher order cor- 
rections to the transition results very tedious, because of 
the large number of coupled flow equations that must be 
analyzed. This is due to the lack of symmetries of this 
model (only parity is preserved), which generates an im- 
portant number of counter-terms that must be taken into 
account in the critical region close to a = 1. In particular 
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we can have higher charges of the form 



M 



h.c.) 



where I > 1, generating next-to- nearest neighbor hop- 
pings and beyond in the low energy theory. However, 
by normal ordering this term can be shown to be irrel- 
evant close to the transition. Also, a renormalization 
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of the potential can be expected, ^ 
v m = V-m hi order to preserve parity symmetry. And fi- 
nally, a renormalization of A which we do not show here. 
The complexity of the problem justifies the application 
of numerical techniques such as the numerical renormal- 
ization group, which deal with the whole effective low- 
energy Hamiltonian. 



B. Variational Ansatz. 

Some deeper insight into the physics of the problem can 
be obtained by using a generalization of the variational 
Ansatz proposed by Silbey and Harris: 27 
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FIG. 1: Phase diagram predicted by the variational calcula- 
tion, for a chain of 5 sites. The model shows a phase transition 
at the critical coupling a c = 1, where the parity symmetry 
is broken. For a < 1 there is a delocalized phase, with t ren 
finite and an effective quadratic potential dependent on the 
coupling strength, that is responsible for the increasing local- 
ization of the particle at the central sites. For a > 1 there 
is a localized phase, with t ren = 0, and a M-fold degenerate 
ground state is predicted. Notice that, in the variational ap- 
proach, a similar phase diagram is predicted for a chain of 6 
sites. 



For the application to our problem we have enlarged the 
original set of variational parameters, fk, in order to in- 
clude the on-site amplitudes of the wave function, c m . 
Without loss of generality, the latter are taken real. The 
energy of the above proposed ground state is: 
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where we have defined t ren = te~? ^ k and g = J\(A+ 
-^) 2 . The minimum condition imposes the following set 
of equations: 

tren(c m +l + C ro _l) + (#(m - TO ) 2 - E G )c m = (7) 



, fk \ 



X 



fk ■ 



((Z 2 ) G (_ + ^) = (8) 



The first set of equations can be seen as those of a particle 
in a chain with renormalized hopping t ren and moving 
in a symmetric parabolic potential v m = g(m — mo) 2 . 
The second equation is a self-consistent condition for the 
parameters fk, once we have determined the lowest Eg 
and the corresponding set of c m from J7])- 



Results 

The solution of these equations gives a phase transition 
for the critical coupling a c = 1, as expected from the pre- 
ceding section. For a < 1 the particle is not localized, 



as t ren is nonzero, but a certain degree of localization at 
the center arises from the existence of the potential term 
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In the limit 



<§; ui c , this factor can be shown to 
t c m+ ic m . As the coupling strength 



be g = 

approaches the critical value, the localization of the par- 
ticle in the central site becomes sharper, as can be seen 
in Fig. [5] for a chain of 5 sites. In this figure the position 
mean square of the particle is plotted as a function of 
the coupling strength. At the critical point the particle 
is completely localized in the center, as a consequence 
of the renormalization to zero of the hopping parameter 
and the appearance of the confining potential. 

For a > 1 the particle gets localized, with t ren — 0, and 
an M-fold degenerate ground state, in which the parity 
symmetry is broken, is predicted. This is, actually, the 
result of solving exactly the Hamiltonian ^ in the lim- 
iting case t — 0, as discussed in Section I. The transition 
line predicted by this approach is vertical, not having any 
dependence on the hopping parameter. A phase diagram 
containing this features is shown in Fig. [TJ 

The variational calculation shows that despite the 
Caldeira - Leggett model includes a counter-term to en- 
sure homogeneous dissipation, the boundary conditions 
are responsible of some non homogeneous effects such as 
the creation of an effective parabolic potential which con- 
fines the particle at the center. However, as we will see in 
the next section, the numerical results suggest that the 
situation close to the critical point is more complicated, 
with new inhomogeneous terms playing an important role 
in the solution. 
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C. Numerical Renormalization Group. 

In order to include higher correlation effects it is nec- 
essary to go beyond the variational solution. A power- 
ful method to study quantum impurity problems is Wil- 
son's numerical renormalization group (NRG) i 28 i 29 Orig- 
inally conceived to deal with fermionic environments, it 
has been recently adapted to handle bosonic bathsj 30 i 31 
Bosonic environments can also be analyzed using the well 
known correspondence between bosons and fcrmions in 
one dimension ) 32 i 33 i 34 i 35 used e.g. in the study of dis- 
sipative gatesi 36 ' 37 However, the fermionic model covers 
only a part of the bosonic one, as the dissipation strength 
is limited to the range < a < 1. Exactly because of 
this fact, we use the bosonic version of the problem. 

When working with bosonic degrees of freedom, we use 
the so called star-NRG 3 -^ method, schematically shown in 
Fig. O As opposed to the conventional "chain" NRG, 
this version allows us to deal with the counter-term in 
the Hamiltonian. This issue arises due to the fact that 
the counter-term, which only involve particle operators, 
is of order of the cutoff. In the chain NRG, the Hamil- 
tonian is transformed in such a way that all the particle 
dependence is included in the first iteration of the algo- 
rithm. The inclusion of the counter-term here requires 
much greater precision in the calculations, as opposed to 
the star-NRG, where is included iteratively. This is no 
longer the case in absence of the counter-term, as shown 
m Ref® There the authors used the chain NRG to 
study dissipative exciton transfer. 

The star-NRG method is based on the introduction of 
a new basis of bosonic states for which the couplings to 
the particle decrease exponentially as A _n , where A is a 
scaling factor between 2 and 3, and n labels the bosonic 
sites. Then, the Hamiltonian is diagonalizcd iteratively, 
adding a new bosonic site in every step. 

The basis is truncated with both an upper cutoff in 
energies, which is progressively reduced in each iteration 
with the scaling factor A, and an upper cutoff in the 
occupation number of the bosonic sites. The first one 
is chosen to have a number of kept states in each iter- 
ation of N s = 100 — 120, while for the second we let 
Ni, = 30 — 40 bosons per site. We check that the proce- 
dure gives well converged results, regardless of the chosen 
truncation parameters. 

We exploit the parity symmetry of the Hamiltonian in 
the code, both for reducing the size of the matrices to be 
diagonalized and for minimizing numerical errors. Let n 
be the parity operator, under which the operators of the 
Hamiltonian transform as: 



UqU^ =-q 



(9) 





FIG. 2: Sketch of the star-NRG Hamiltonian used in the work, 
a) Every bosonic site a n couples to the particle. In each it- 
eration a new site is added and the resulting Hamiltonian 
is diagonalized, giving the energy spectrum, b) The particle 
Hamiltonian, coupled to the first bosonic site. Notice that the 
structure of the couplings is the same for the rest of bosons. 
See the text for more details on the calculations. 



the basis of eigenstates of parity, c„ 
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For M odd the central site is always guaranteed to be 
a parity eigenstate. The states used to diagonalize the 
NRG-Hamiltonian at zero iteration are then \m,p; n&; P), 
with total parity P = p(—l) nb being a good quantum 
number. The total Hamiltonian then splits into two sep- 
arated sectors of well defined total parity, reducing by 
two the size of the matrices to be diagonalized. The same 
can be shown to be true at iteration N + 1, where the 
total parity states \r,p; rib; P)n+i are constructed adding 
a new bosonic site to the eigenstates at iteration N with 
parity p. The matrix elements of the NRG-Hamiltonian 
verify, in this basis: 



(r',p';ri b ;P'\HN+i\r,p;n b ;P) cx 5 



Bosonic states have a well defined parity, n(6j ; )™''|0} = 
(-i) nb (bl) n »\Q). The particle states must be rotated to 



>p,p' 



The output of the NRG procedure are the flows of the 
lowest lying energy states as the cutoff is reduced itera- 
tively. At some point the flows are expected to converge 
to stable (low energy) fixed points. The effective Hamil- 
tonian can be reconstructed analyzing the evolution of 
those flows, as well as the evolution of other observables 
of the system. Here we will use the evolution of the aver- 
aged position of the particle, (q)N, and its mean squared 
deviation (q 2 )N, evaluated in the ground state. Those 
flows are enough to characterize the different phases of 
the system. 
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FIG. 3: Representative flows of the Numerical Renormal- 
ization Group transformations carried out in this work, for 
M = 5 sites. The horizontal axis is the iteration number, and 
the graphs are scaled energy levels, a) Flow towards a non 
degenerate ground state (a — 0.8, t = 0.01) The inset shows 
the flow of the mean squared position of the particle, which 
becomes localized around the center, b) Flow towards a de- 
generate ground state (a = 1.2, t = 0.01). The lower inset 
gives details of the way in which the lowest energy levels flow 
to the fixed point. The top inset shows the localization of the 
particle beyond the center, as indicated by its mean squared 
position. Notice that, in this regime, two energy scales are 
playing a part in the flow, delimitating an intermediate regime 
which is dominated by an unstable fixed point. 



Results 

In this paper we analyze chains with a number of sites 
ranging between 2 and 6. For two sites the NRG repro- 
duces the phase diagram of the dissipative two level sys- 
tem, as was shown by Bulla et al^i. For larger chains and 
small dissipation, the results are in quantitative agree- 
ment with the variational solution (see below), predict- 
ing a delocalized phase with renormalized hopping and a 
renormalized potential which tends to localized the par- 
ticle at the center as the coupling strength is enlarged. 
As far as the phase transition is concerned, the case of 
AI = 3,4 does not deviate too much from the dissipa- 
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FIG. 4: Flows of the NRG for M = 6 sites, a) Flows in 
the first region of the localized phase (a = 1.2, t = 0.01), 
where the particle is confined at the central sites, b) Flow in 
the second region of the localized phase (a = 1.4, t = 0.01). 
Here the particle is confined in the next to the center sites, 
suggesting the formation of a double well effective potential 
in the array. In both cases the lower insets show in detail the 
flow of the lowest energy levels. Notice that, as happened in 
the chain of five sites, two energy scales apparently delimit the 
onset of an unstable fixed point, corresponding to an effective 
chain of M — 4 sites. This is actually well conformed by the 
value of the mean squared position of the particle (inset) 



tive two level system. In both cases there is a phase 
transition in which t ren = 0, but for M = 3 the parity 
symmetry is not broken, because the particle is localized 
at the center. For M = 4 the phase transition is that 
of the dissipative two level system, the edge sites being 
decoupled in energy from the central ones. This is all in 
contrast with the variational solution, where a transition 
to an M degenerate state is predicted. 

Of more interest are the cases of M = 5,6. Here a 
new behavior is observed, which should be representa- 
tive of the one expected for larger chains. The energy 
flows for small and large dissipation are shown in Fig. [3] 
and Fig. 2] Again, weak dissipation induces some local- 
ization at the center of the array of the particle density, 
as can be seen in the inset of the figures, where (q 2 )at 
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FIG. 5: Mean-squared position of the particle as a function of 
the coupling strength in the delocalized phase of the model, 
for M = 5, as predicted by the variational calculation and 
the NRG. Both approaches show an increasing localization of 
the particle at the center of the chain as the coupling strength 
gets larger. This effect arises due to the renormalization of the 
hopping parameter and the emergence of an effective confin- 
ing potential. In the variational calculation the localization is 
stronger, suggesting that higher order corrections to the effec- 
tive potential are playing an important role in the numerical 
calculation. 



is computed. As mentioned above, in this regime the 
results agree qualitatively with the variational solution, 
as shown in Fig. where the mean squared position of 
the particle is calculated in both approximations for a 
chain of 5 sites (a similar plot can be obtained for 6 sites, 
the main difference being that the mean squared position 
tends to a finite value for a — * 1, having two states in 
the center instead of one) . 

The differences between the localized phase predicted 
by the variational calculation and the NRG are even 
sharper for those longer chains. Above a critical strength 
coupling a c of order one, we find a doubly degenerate 
state, for odd and even number of sites, in which the 
parity symmetry is broken. For M — 5 the particle lo- 
calizes next to the center, while for M = 6 it does initially 
in the central sites, and for larger dissipation in the next 
to the center ones. This result follows from analyzing the 
degeneracy of the ground state, extracted from the en- 
ergy flows, as well as the evolution of the mean squared 
position operator. The converged values of the latter can 
be used to make an ansatz of the sort of ground state 
density matrix to which the flow converges. 

In Tablefl] the zoo of stable fixed points of the model is 
shown, for chains of M = 5,6 sites. In the case of M = 6 
there is an extra fixed point in the localized regime, cor- 
responding to a situation in which the particle finds more 
favorably to get localized in the sites next to the edges 
than next to the central sites. This second transition also 
occurs for a critical value of the coupling strength, but 
there is neither a symmetry breaking nor a change in the 
degeneracy of the ground state. Thus, the information 
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TABLE I: Stable fixed points of the NRG for chains of M = 
5, 6 sites. The fixed points are characterized by the ground 
state degeneracy, GSdeg and the position mean squared value, 
(q 2 ) nrg, whose values can be obtained with the NRG (in the 
localized phase there is a single value for every a, while in the 
delocalized one the value depends on the coupling strength, as 
shown in Fig. 5 for a chain of 5 sites) . Those are used to pro- 
pose an ansatz for the ground state density matrix, p, which 
fits it correctly. The states \i) represent a particle sitting at 
site i. \t ren ,a) is the ground state of a free tight-binding 
chain with hopping t ren and a parabolic potential, dependent 
on the strength coupling a. The latter is the output of the 
variational calculation for the delocalized phase, which fits 
quite well the numerical data. 

provided by NRG is not enough to fully characterize the 
nature of this phase transition. 

This is not the only limitation of the numerical 
method. It also does not allow us to study large val- 
ues of the dissipation, a ^> 1, as the occupancies of the 
bosonic states become also high. The question of whether 
other phase transitions can be ruled out for high values 
of a remains open, and deserves a separate study with 
different techniques. 
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FIG. 6: Plot of the dependence of the crossover scale T* on the 
distance to the critical coupling strength. Here T* oc A _JV , 
with JV* chosen as the iteration for which the first excited 
level verifies A N En*,i = 0.03. The figure shows the results 
for two different hopping parameters. In both cases there is a 
good agreement with an exponential decay of T* as a function 
of the distance from the critical coupling, log T* oc l/(a c — a) 

From the energy flows some extra information can be 
extracted. In the delocalized phase, a single energy scale 
seems to be playing a part in the evolution from high 
energy to low energy behavior. Actually, the flow in this 
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phase is similar to that in the dissipative two level system, 
and in the same way we can define a crossover scale T* oc 
A~ N from the iteration N* at which the flow changes 
from its initial behavior to the low energy regime. As 
shown in Fig. [6] T* tends to zero exponentially as the 
coupling strength approaches the critical value, log T* oc 
\/{a. c — a). Hence, our results suggest that the transition 
is continuous, being consistent with the existence of a KT 
transition^ 

The flows in the localized phase show a different be- 
havior. Here, two different energy scales appear in the 
course of the flow, revealing an unstable fixed point in 
an intermediate regime. Those scales are defined now 
from the iterations iV-j* and N£ at which the energy lev- 
els decouple from the low energy sector, giving rise to two 



crossover temperatures Ti oc A 



-NT 



1,2. From the 



mean squared position flows it can be deduced that the 
upper energy scale marks the decoupling of the sites lo- 
cated at the edges, as the values of this operator are well 
fitted to the expected ones in free tight-binding chains 
with finite hoppings but sites in the edges supressed (ef- 
fectively reducing the chain to one with two less sites). In 
this way, the intermediate fixed point would correspond 
to an effective cluster of three sites in the M = 5 case, 
and four sites in the M = 6 one, with a renormalized 
hopping parameter t ren . The lower energy scale corre- 
sponds to the onset of the phase transition, for here the 
parity symmetry is broken and only two sites remain in 
the low energy regime. 
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IV. PHASE DIAGRAM 

A phase diagram of the model, including all the fea- 
tures discussed above, is presented in Fig. [7] As in the 
widely studied dissipative two level system, there is a 
phase transition between a delocalized regime to a local- 
ized one. In the delocalized phase the effect of the bath is 
that of reducing the effective hopping and of generating 
a renormalized potential which makes the density of the 
particle higher close to the center. In the localized phase 
the parity symmetry is broken and in both cases, odd 
and even, the particle localizes in one of two degenerate 
sites. This transition is continuous, and the numerical 
results are consistent with a transition of the KT type, 
as in similar dissipative systems. 

The variational calculation shows how a parabolic ef- 
fective potential emerges from the coupling to the bath, 
being responsible of some localization of the particle at 
the center. However, the numerical results suggest the 
existence of a more complicated renormalized potential 
which would explain both the almost complete localiza- 
tion of the particle as a — > 1 from below, and the in- 
homogeneous degenerate ground state in the localized 
phase. A simple guess which works well qualitatively is a 
effective potential in the form V m = (<7 /2)(to — mo) 2 + 
(<7i/4!)(m — mo) 4 - If ,9o oc ip-c ~ a )i this Ansatz has 
a minimum at m — too in the delocalized phase, and 



FIG. 7: Phase diagram of the model, eq.([T|, for M = 5 (top) 
and M — 6 (bottom) sites, deduced from the NRG flow. The 
continuous lines are linear fittings to the numerical data. The 
insets show representative density distributions of the particle 
in each phase. In both cases there is a phase transition to 
a localized phase in which the parity symmetry is broken. 
The even case shows also a second transition in the localized 
phase, where the particle is confined beyond the central sites, 
resembling the kind of localization observed in the odd case. 



at to = too ± w gi 9 ° in the localized one, explaining 

not only the doubly degenerate state, but also that in 
the case M = 6 there is a second transition to a phase 
in which the particle localizes in the next to the center 
sites. From the point of view of renormalization theory, 
such an ansatz is reasonable as higher corrections to the 
potential should be more irrelevant. 

In the region of the phase diagram which we have an- 
alyzed (a < 2), we have not found any further crossover 
to a region where the particle is confined to the edges. 
This could be explained by the role played by the in- 
termediate unstable fixed point in the localized phase. 
By studying the case of three and four sites, the only 
way to get a phase transition in which the particle lo- 
calizes at the edges is by starting with slightly lower site 
energies here as compared to the center. Hence, the de- 
coupling of the edges would be necessary to give rise to 
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such a renormalization of the on-site energies, favoring 
the phase transition to a more stable regime. 

For larger chains we expect a similar picture to work, 
and new features should not appear as far as the con- 
tinuous limit is not reached. The phase diagram should 
show a transition to a localized phase for a critical value 
a c also around the unity. Close to this phase transition, 
the renormalized potential gets quartic corrections and 
the particle becomes localized in the resulting double- 
well profile. As the minimum of this effective potential 
depends on the dissipation strength, there should be sev- 
eral crossovers to regions in which the particle is localized 
at points increasingly farther from the center. However, 
the particle is never localized at the edges, since their 
decoupling seems to be crucial to the realization of this 
inhomogenous transition. Thus, in the localized phase, 
two or more energy scales are expected to play a role, 
depending on the number of energy levels that are de- 
coupled from the low-energy sector until the stable fixed 
point is reached. In the continuous model, which corre- 
sponds to the case of an infinite number of sites in the 
array, we expect only a single phase transition, at a c 
around the unity, to a phase where the particle is local- 
ized in a double well potential profile whose minimum 
depends on the coupling strength. 

V. CONCLUSIONS. 

We have analyzed the simplest extension of the widely 
studied dissipative two level system, which interpolates 
between this model and the also extensively analyzed dis- 
sipative quantum particle in a periodic potential. The 
main results are obtained with a numerical renormaliza- 
tion group technique specifically adapted to bosonic co- 
ordinates. For small and intermediate values of the dissi- 
pation parameter, a, our results provide a well controlled 
approximation to the ground state of the system. 

When dissipation is weak, we find that the density dis- 
tribution of the particle becomes narrower, and localized 
around the center of the array. This result is consis- 
tent with the expectation that the environment acts as 
a measurement apparatus on the position of the parti- 
cle, leading to a more localized distribution. The initial 



degeneracy between the sites of the array is broken by 
the combination of the abrupt boundary conditions and 
the dissipation. As a result, an effective local potential 
is induced, with a minimum at the center of the array, 
leading to the confinement of the particle near the center. 

For values of the dimensionless dissipation strength 
a ^ 1, we find a transition to a situation with a doubly 
degenerate ground state, in which the parity symmetry 
is broken. In this case, the combined effect of boundary 
conditions and dissipation leads to the formation of an 
effective double well potential irrespective of the number 
of sites in the array. 

The formation of an inhomogeneous effective potential 
from a spatially homogeneous coupling to a dissipative 
bath should be a generic feature in similar models. A 
more remarkable result is that, in some cases, this effec- 
tive potential is non-monotonic, so that simple confine- 
ment geometries can lead to complicated patterns in the 
localized regime for high values of the dissipation. An 
open question is to what extent these results may de- 
pend on the particular choice of particle-bath coupling. 
It has been argued in Ref. [39| that the introduction of 
a counter-term is insufficient to introduce full transla- 
tional invariance in some dynamic contexts, such as that 
of a suddenly introduced coupling. An extension of the 
present equilibrium study to models of truly translation- 
ally invariant dissipation 39 could shed some light on this 
issue. 
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